超临界压力CO2在水平圆管内流动传热数值分析

采用SST k-w低雷诺数湍流模型对加热条件下超临界压力CO2在内径di=22.14 mm,加热长度Lh=2440 mm水平圆管内三维稳态流动与传热特性进行了数值计算.通过超临界CO2在水平圆管内的流动传热实验数据验证了数值模型的可靠性和准确性.首先,研究了超临界压力CO2在水平圆管内的流动传热特点,基于超临界CO2在类临界温度Tpc处发生类液-类气“相变”的假设,揭示了水平圆管顶母线和底母线区域不同的流动传热行为.然后,分析了热流密度qw和质量流速G对水平圆管内超临界压力CO2流动换热的影响,通过获取流体域内的物性分布、速度分布和湍流分布等详细信息,重点解释了不同热流密度qw和质量流速G下顶母线内壁温度Tw,i分布产生差异的传热机理,分析结果确定了类气膜厚度d、类气膜性质、轴向速度u和湍动能k是影响顶母线壁温分布差异的主要因素.研究结果可以为超临界压力CO2换热装置的优化设计和安全运行提供理论指导.     

水化硅酸钙纳米压痕分子动力学模拟方法优化

针对水化硅酸钙纳米压痕模型忽视了压头与基底之间相互作用的问题，由尺寸差异引起的金刚石压头难以计算的问题，以及Wittmann模型无法得到实际接触面积的问题，提出了新的模型与计算方法．结合分子动力学方法，采用金刚石压头-Wittmann模型基底的组合方式构建无定形态水化硅酸钙纳米压痕试验模型．在建模阶段，考虑到压头模型与基底模型粒子间尺寸差异，提出了等比例替换模型，通过公式推导并就不同尺寸模拟结果验证了等比例替换模型的可行性．在计算阶段，提出了局部前处理的弛豫方法进行模拟．确定最大荷载位置处的接触面积为546 nm2，进而求出水化硅酸钙模型硬度H为0.84 GPa、折合模量Er为30.52 GPa．并通过纳米压痕试验，验证了模拟结果的准确性，证明了模型的科学性，对今后水化硅酸钙(C-S-H)纳米层面的模拟具有重要借鉴意义．     

染料敏化太阳能电池载流子传输的数值模拟

由于在染料敏化太阳能电池(dye-sensitized solar cell, DSSC)中存在染料弛豫、半导体薄膜中电子与氧化态染料分子发生反应和电子在电解质中与氧化态离子复合等不利反应,利用一个更完善的DSSC载流子传输模型对电池的光电性能进行模拟就显得非常重要。为此,本文基于由多重俘获理论建立的DSSC中的包括电子、染料阳离子、碘化物和三碘化物在内的载流子传输模型,数值模拟得到了不同TiO₂薄膜厚度、不同入射光强度与不同染料分子吸收系数下DSSC的J-V曲线。结果表明,随着TiO₂薄膜厚度的增加,太阳能电池的短路电流密度增大,开路电压减小,光电转换效率先增大后减小。当DSSC的TiO₂薄膜厚度为20 μm时,光电转换效率达到最大值7.41%,同时光电转换效率随入射光强度与染料分子吸收系数的增大均有一定程度提高,其中在吸收系数为4 500 cm^-1时,光电转换效率为6.73%。以上结果可以为改进DSSC的光电性能提供理论指导。     

Numerical simulation of transpiration cooling through porous material

Transpiration cooling using ceramic matrix composite materials is an innovative concept for cooling rocket thrust chambers. The coolant (air) is driven through the porous material by a pressure difference between the coolant reservoir and the turbulent hot gas flow. The effectiveness of such cooling strategies relies on a proper choice of the involved process parameters such as injection pressure, blowing ratios, and material structure parameters, to name only a few. In view of the limited experimental access to the subtle processes occurring at the interface between hot gas flow and porous medium, reliable and accurate simulations become an increasingly important design tool. In order to facilitate such numerical simulations for a carbon/carbon material mounted in the side wall of a hot gas channel that are able to capture a spatially varying interplay between the hot gas flow and the coolant at the interface, we formulate a model for the porous medium flow of Darcy–Forchheimer type. A finite‐element solver for the corresponding porous medium flow is presented and coupled with a finite‐volume solver for the compressible Reynolds‐averaged Navier–Stokes equations. The two‐dimensional and three‐dimensional results at Mach number Ma = 0.5 and hot gas temperature T_HG=540 K for different blowing ratios are compared with experimental data. Copyright © 2014 John Wiley & Sons, Ltd.     

Blended-fuel based EDC combustion model and its application in heptane-ethanol fire simulation

The blended-fuel based eddy-dissipation-concept combustion model was newly developed in the FireFOAM framework, and applied to simulate 30 cm×30 cm heptane-ethanol pool fire. Comparison was made of fire height, centerline temperature against experimental measurements, which shows that they match very well with each other. However, further studies are needed to examine the validation of this model in fire simulations with various scales.     

In Situ Detection of Particle Aggregation on Electrode Surfaces

Partially blocked electrodes (PBEs) are important; many applications use non‐conductive nanoparticles (NPs) to introduce new electrode functionalities. As aggregation is a problem in NP immobilization, developing an in situ method to detect aggregation is vital to characterise such modified electrodes. We present chronoamperometry as a method for detection of NP surface aggregation and semi‐quantitative sizing of the formed aggregates, based on the diffusion limited current measured at PBEs as compared with the values calculated numerically for different blocking feature sizes. In contrast to voltammetry, no approximations on electrode kinetics are needed, making chronoamperometry a more general and reliable method. Sizing is shown for two modification methods. Upon drop casting, significant aggregation is observed, while it is minimized in electrophoretic NP deposition. The aggregate sizes determined are in semi‐quantitative agreement with ex situ microscopic analysis of the PBEs.     

Error Estimates for a Class of Elliptic Optimal Control Problems

In this article, functional type a posteriori error estimates are presented for a certain class of optimal control problems with elliptic partial differential equation constraints. It is assumed that in the cost functional the state is measured in terms of the energy norm generated by the state equation. The functional a posteriori error estimates developed by Repin in the late 1990s are applied to estimate the cost function value from both sides without requiring the exact solution of the state equation. Moreover, a lower bound for the minimal cost functional value is derived. A meaningful error quantity coinciding with the gap between the cost functional values of an arbitrary admissible control and the optimal control is introduced. This error quantity can be estimated from both sides using the estimates for the cost functional value. The theoretical results are confirmed by numerical tests.     

A modified pseudospectral method for solving trajectory optimization problems with singular arc

This paper presents a direct method based on Legendre–Radau pseudospectral method for efficient and accurate solution of a class of singular optimal control problems. In this scheme, based on a priori knowledge of control, the problem is transformed to a multidomain formulation, in which the switching points appear as unknown parameters. Then, by utilizing Legendre‐Radau pseudospectral method, a nonlinear programming problem is derived which can be solved by the well‐developed parameter optimization algorithms. The main advantages of the present method are its superior accuracy and ability to capture the switching times. Accuracy and performance of the proposed method are examined by means of some numerical experiments. Copyright © 2016 John Wiley & Sons, Ltd.     

Effects of awareness program and delay in the epidemic outbreak

In the present paper, an epidemic model has been proposed and analyzed to investigate the impact of awareness program and reporting delay in the epidemic outbreak. Awareness programs induce behavioral changes within the population, and divide the susceptible class into two subclasses, aware susceptible and unaware susceptible. The existence and the stability criteria of the equilibrium points are obtained in terms of the basic reproduction number. Considering time delay as the bifurcating parameter, the Hopf bifurcation analysis has been performed around the endemic equilibrium. The direction of the Hopf bifurcation and the stability of the bifurcating periodic solutions are investigated by using the normal form theory and central manifold theorem. To verify the analytical results, comprehensive numerical simulations are carried out. Copyright © 2016 John Wiley & Sons, Ltd.